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ABSTRACT 

Cosmic ray particles are more rapidly accelerated in oblique shocks, with the magnetic 
field inclined with respect to the shock normal direction, than in parallel shocks, as a 
result of mirror reflection at the shock surface and slower diffusion in the shock nor- 
mal direction. We investigate quantitatively how these effects contribute to reducing 
the acceleration time over the whole range of magnetic field inclinations, ft is shown 
that, for quasi-perpendicular inclination, the mirror effect plays a remarkable role in 
reducing the acceleration time; whereas, at relatively small inclination, the anisotropic 
diffusion effect is dominant in reducing that time. These results are important for a 
detailed understanding of the mechanism of particle acceleration by an oblique shock 
in space and hcliospheric plasmas. 

Key words: acceleration of particles — diffusion — magnetic fields — shock waves 
— cosmic rays. 



1 INTRODUCTION 

Diffusive shock acceleration (DSA) is one of the most 
favourable mechanisms to reproduce the energy spectrum of 
observed cosmic rays (for a review, see Drury 1983). Shock 
waves accompanied by the magnetic fields with small fluctu- 
ations (Alfven waves) are well established in space plasmas 
and operate as a powerful accelerator of charged particles 
jBlandford fc Eichlerlll987l) . In earlier works, substantial ef- 
forts were devoted to studies of acceleration by the sim- 
ple 'parallel shocks' in which both the magnetic field and 
the direction of pl asma flow are parallel to the direction of 
shock propagation jAxford. Leer fc Skadronll97ctlBeljl973 
iBlandford fc Ostrikerll978lb At the parallel shock fronts, all 
upstream particles are advected into the downstream region, 
gyrating around the mean magnetic field. According to the 
conventional resonant scattering theory, a particle with its 
gyroradius comparable to the turbulent Alfven wavelength 
is resonantly scattered and migrates back upstream after 
many times of small-angle scattering. Cosmic ray particles 
acquire their energies in the process of transmission back and 
forth across the shock. Since the particles gain only a small 
amount of kinetic energy in each traversal of the shock, the 
acceleration efficiency depends largely on the rate of shock 
crossing of the particles. 

Relating to particle acceleration in the solar wind and 
the Earth's bow shock, the more general 'oblique shocks', 
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across magnetic fi eld lines, have been studied llToptvghinl 
ll980tlDeckeilll98ot . Some researchers have argued that the 
acceleration efficiency is enhanced in oblique shocks, com- 
pared with that in parallel shocks <IJokipiilll987t lOstrowskil 
1988). If the particle momenta along the magnetic field, 
Pll , are larger than a critical value, the particles gain en- 
ergy solely via successive transmissions through the shock 
front, just as in the case of parallel shocks. However, when 
the value of p\\ is smaller than the critical value, upstream 
particles cannot penetrate into the downstream region with 
stronger magnetic field and are reflected back into the up- 
stream region with weaker field, having their pitch angles 
reversed. On the basis of the conservation of magnetic mo- 
ment, the turnover point of the pitch angle is determined 
by the ratio of the upstream/downstream magnetic field 
strength. This 'mirror reflection' results in significant re- 
duction of acceleration time. Quite efficient acceleration is 
expected, in particular, for quasi-perpendicular shocks with 
larger inclination of magnetic field lines. 

It was also pointed out that the acceleration time for 
oblique shocks could be reduced owing to the anisotropy of 
the particle diffusion coefficient (|Jokipiilll987tb The effective 
diffusion coefficient involved in the DSA time-scale can be 
represented by the tensor component normal to the shock 
surface, K n , which is decomposed into components parallel 
(/t||) and perpendicular (ki) to the magnetic field (see also 
Section 2.1). In the special case of parallel shocks, the effec- 
tive diffusion coefficient reduces to K n = Ka . In the ordinary 
case of fey 2> reflecting one-dimensional (ID) properties 
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in th e conventional m agnetohydrodynamic (MHD) descrip- 
tion llOstrowskilll988l) . the value of k u decreases as the mag- 
netic field inclination increases, and in perpendicular shocks, 
K n — k±. Within the DSA framework, the smaller value of 
«n leads to a shorter acceleration time, and thereby a higher 
acceleration efficiency. 

Summarizing the above discussions, in oblique shocks 
there exist two possible effects contributing to the reduction 
of acceleration time: mirror reflection at the shock surface, 
and diffusion anisotropy. In previous works, these contribu- 
tions were comprehensively treated with the 'oblique effect'. 
In the present paper, we investigate the contribution of these 
effects separately and reveal which of these two effects con- 
tributes more dominantly to the acceleration time over the 
whole range of field inclination. For this purpose, we derive 
the expression for the acceleration time with and without 
the effects of mirror reflection. Presuming that no particles 
are reflected at the oblique shock, we estimate the acceler- 
ation time without mirror effects, though still including the 
effects of diffusion anisotropy. To our knowledge, so far there 
has been no publication that quantitatively reveals the ef- 
fects of anisotropy on the reduction of acceleration time for 
oblique shocks. Here we demonstrate that mirror reflection 
makes a dominant contribution to the rapid acceleration for 
highly inclined magnetic fields, whereas anisotropic diffusion 
becomes effective at smaller inclination angles. 



2 TIME-SCALES OF OBLIQUE SHOCK 
ACCELERATION 

In the following description, we use a unit system in which 
the speed of light is c = 1. Upstream and downstream quan- 
tities are denoted by the indices i=l and 2, respectively. The 
subscripts || and _!_ refer to the mean magnetic field, and sub- 
script n refers to the shock normal. The subscripts a indicate 
the processes of particle interaction with the shock: a =r for 
reflection and a —12 (21) for transmission of the particles 
from region I (2) to region 2 (I). For non-relativistic flow 
velocities, calculations related to particle energy gain are 
performed up to first order in Vi, where Vi = Ui/ cos 8i, and 
Ui and 9i are the flow velocity in the shock rest frame and 
the magnetic field inclination to the shock normal, respec- 
tively. 

2.1 Effective diffusion coefficient 

To evaluate the DSA time-scale of cosmic rays, we need their 
spatial diffusion coefficient. According to the tensor trans- 
formation, the effective diffusion coefficient for th e shock 
normal direction can be expressed as dJokipiilll987l) 

Kni = COS 2 9i + Kl, SU1 2 6i. (1) 

Without going into the details of the sc attering process, here 
we make use of the empirical scalings (lOstrowskilll 988T) 

/tii = kbx and k± = kb/x. (2) 

Here, the scaling factor I) reflects the energy density 
of the mean to fluctuating magnetic fields, (B/SB) 2 , and 
kb = r s v/3 is the Bohm diffusion coefficient, where r g and 
ii are the gyroradius and speed of a test particle, respectively. 



Making the assumption that the magnetic moment of 
the particle is approximately conserved during the inter- 
action with the discontinuous magnetic field at the shock 
front, we consider the situation in which small fluctuations 
are superimposed on the mean regular field, i.e. x ^> 1, 
in the vicinity of the shock (free crossing limit). On the 
other hand, in the case of x ~ 1 (diffusive limit), the 
fluctuations significantly affect the coherent gyromotion of 
the particle, and precise numerical simulations are required 
JPecker fe Vlahosl 1 1986ft . For x 1, we have the relation 
of k±(= k\\/x 2 ) <C Kn. This means that, for large angles 
6i > tan -1 x, the term involving k± becomes dominant on 
the right-hand side of equation (1). Thus, one finds that 
larger values of x and Qi are likely to lead to a shorter ac- 
celeration time. Here it is noted that the acceleration time 
cannot be infinitesimally small, because of the limits of K n 
and 9i (see Section 2.2). 

The parallel/perpendicular diffusion of cosmic rays 
in turbulent plasmas is a long-standing issue and is not 
yet completely resolved. Analytic models are consistent 
with neither observations (iMazur et al. | 2000; Zhan g et alJ 
2003) nor numerical simulat ions (Giac alone fc Jokipiilll999l : 
iMace. Matthaeus fc Bieb er 2000) over a wide energy range. 
This situation has changed since the dev elopment of non- 
linear guiding centre theory (NLGC) by iMatthaeus et all 
(2003). Recent calculation of the NLGC diffusion coefficient 
suggests reasonab le values of x 2 > 10 3 for the typical solar 
wind parameters JZank et al.ll2004ft . which can be accomo- 
dated with the free crossing limit [x 1). In many astro- 
physical environments, the perpendicular diffusion of cosmic 
rays still remains uncertain and is commonly inferred from 
fluctuations of the turbulent magnetic field. We will discuss 
this point in the conclusions. 

As long as the values of x can be regarded as con- 
stants, the present consequences concerning the accelera- 
tion time are generic, apparently independent of the type 
of turbulence. However, for an application, when evaluating 
the maximum possible energy of a particle by equating the 
acceleration time with the competitive cooling time-scales, 
one should pay attention to the fact that, in general, x de- 
pends on the ratio of r g to the correlation length of tur- 
bulence, deviating from the simplistic scaling x ~ (B/SB) 2 
llHonda fc Hondall2004l) . Here, note that, in the Bohm limit 
with v = 1, the above dependence is found to be weak, so as 
to appear as logarithmic. For the present purpose, below, we 
assume i to be a constant much larger than unity, reflecting 
weak turbulence, and thereby larger anisotropy of diffusion. 



2.2 Particle acceleration including mirror effects 

Follow ing the kinematic approach developed by lOstrowskil 
(1988), for instruction we outline the derivation of the ac- 
celeration time for an oblique shock including the effects 
of both mirror reflection and diffusion anisotropy. For con- 
venience, the energy and momentum of a test particle are 
transformed from the upstream/down stream rest frames to 
the de Hoffmann- T eller (HT) frame dde Hoffmann fc Teller! 
ll950UHudsoni ri965). Since the particle kinetic energy is in- 
variant in the HT frame, where electric fields vanish, one 
can easily estimate the energy gain of the particle during 
the interaction with the shock, presuming the conservation 
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of magnetic moment of the particle. Finally, all variables are 
transformed back to the original rest frames. 

In the HT frame, if the cosine of particle pitch angle, 
H, is in the range < n < (io, the upstream particles are 
reflected at the shock surface. Here, 



Mo 



•fli/Ba) 



1/2 



gives the critical value, where Bi < B2 for the fast-mode 
oblique shock. In this case we have the following ratio of 
particle kinetic energies: 

^=7i 2 (l + V 1 2 + 2^ lM ), 

where E and E r are the particle energies (in the region 1 
rest frame) before and after reflection, respectively, and 71 = 
(1 — V] 2 ) -1 / 2 . If /lo < /1 ^ 1, the particles are transmitted 
to region 2. In this case, we have 



-Ei 2 
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7172 < 1 + Vifi - V 2 



(1 + Virf 



7i 2 (l~M 2 ) 
(1 - 14) 
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where E and E\ 2 are the particle energies (in the region 
1 and 2 rest frames) before and after transmission, respec- 
tively, and 72 = (1 — Vsj 2 ) -1 . For the transmission of par- 
ticles from region 2 to region 1, 
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where E21 the particle energy (in the region 1 rest frame) 
after transmission. 

The mean acceleration time is denned as t\ = At/d, 
where At is the cycle time and d the mean energy gain per 
interaction. Ignoring particle escape, the cycle time can be 
written in the form 



At = tlP r + (tl +t 2 )Pl2, 



(3) 



where t;[= 2« n i / (v n iUi)] denotes the mean residence time in 
region i and Pa is the probability for process a. Note that, 
for K|| ^> k±, the mean normal velocity can be estimated as 
v n i — «|| yj K n i/K\\i ~ «|| cos6*i, where v\\ = v/2 ~ c/2. 

The probabilities of reflection and transmission are ex- 
pressed as P r = Si /(St + S12) and Pi 2 = Sia/(& + S12), 
respectively. Here S a denotes the normal component of par- 
ticle flux flowing into the shock surface, which is calculated 
by using the normal velocity of the guiding centre drift of 
the particles, that is, 

V le ldn, 

where 

Kd = (mcos6»! + - MCOS^lf/l) 

for a =12 and r (in the region 1 rest frame) and 

V te i = (fiCOs6 2 ~ U 2 )/(l + VCOS62U2) 

for a =21 (in the region 2 rest frame). Carrying out the 
integrations, one approximately obtains P r ~ /x 2 , and P12 — 
1 — /i 2 . Using equations (1) and (2), the cycle time (equation 
3) can then be expressed as 



At 



2kbx 
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I 2 a , sin" 0i , r cos'ffi + ^ sm^i I 
X I cos ei + _^ + (cos2gi+r2sin2gi)3/2 | (4 ) 

where r = U1/U2 is the shock compression ratio. 

The mean energy gain of the upstream particle is de- 
noted as 



d = d r P r + (di2+d 2 i)Pi2, 
where 

d a = / V rcl {E a I E - 1) d^/S a 



(5) 



defines the fractional energy gain in the process a. To an 
approximation up to first order in Vi, then, the resultant 
expressions are 



r oVi 



and 

dl2 = d2i 



■ t [V 1 (l-^)/(l-^ )-V2]. 



Substituting these results into equation (5) gives 
d~ 1^-^(1 -Mo)]- 

As a result, we arrive at the following expression for the 
mean acceleration time (Kobayakawa, Honda & Samura 
2002): 



t A 



2, , sin 2 6>i r(cos 
cos 0i H h 



+ 



(cos 2 



^3/2 



(6) 



Here we have replaced all the downstream quantities with 
upstream ones. Note that equation (6) is valid for the free 
crossing limit x 3> 1. 

In the allowed range of magnetic field inclination an- 
gles, @i ^ cos -1 U\ (de Hoffmann & Teller 1950), the value 
of tA(#i 7^ 0°) is smaller than that of tA (0i =0°). Relating 
to the reduction of t\ , we note that n n involved in the ac- 
celeration time (equation 6) can take a value in the range 
> Uif /\V f\, where / is the phase space density of particles. 
This inequality just reflects the condition that the diffusion 
velocity of particles must be larger than the shock speed. 
Thus, the requirement that the gyroradius cannot exceed 
the characteristic length of the density gra dient (//|Vf |), 
recasts the above condition into K n > r g Ui ljokipiilll987fi . 



2.3 Particle acceleration without mirror effects 

Following the procedure explained above, we derive the ex- 
pression for the mean acceleration time, excluding mirror ef- 
fects. Assuming that all upstream particles are transmitted 
downstream through the shock front (no reflection), thereby, 
setting P r = and P12 = 1 in equation (3), reduces the ex- 
pression of the cycle time to 



At' = ti +t 2 . 

Equation (7) can be explicitly written as 
2kbx 



(7) 
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2 sin 0i 
cos 0i H ^ — 



+ (r 2 /a; 2 ) sin 2 



(8) 



for x 2> 1. Note that the difference from equation (4) is only 
the denominator of the third term in the curly brackets. 

Similarly, the mean energy gain per interaction is de- 
noted as 



d! = d r2 + d 2 



(9) 



Recalling that all particles having pitch angle of < fj. ^ 1 
are forced to be transmitted to region 2, we have 

d'~l(f/i-f7 2 ). 

In contrast to d(oc l/cos0i), the expression for d' is in- 
dependent of the field inclination, and appears to be iden- 
tical to that for parallel shocks. Using equations (8) and 
(9), the acceleration time without mirror effects, defined by 
4 = At' /d', is found to be represented in the following 
form: 



t'A = 



U 2 (r - 1) cos 0i 



2 sin fi 
cos 0i H — 



+ (r 2 /x 2 )i 



10) 



In comparison with equation (6), the value of equation 
(10) is boosted as a result of the factor (cos0i) _1 and the 
smaller denominator of the third term in the curly brackets. 
The latter comes directly from exclusion of mirror effects. 
The factor (cos0i)~ reflects the anisotropy of particle ve- 
locity, involved in the mean residence time tj. Although in 
the evaluation of tA this factor was canceled out by the same 
factor from the expression for d, in the present case it is not 
cancelled, because of the independence of 0i in d' . Note the 
relation 

4(<?i = 0°) = U(0i = 0°) = (3K B x/U?)[r(r + l)/(r - 1)], 

which coincides with the acceleration time for a parallel 
shock. 



3 NUMERICAL RESULTS 

Below, we assume a monatomic, non-relativistic gas with 
specific heat ratio of 5/3, whereby 1 < r ^ 4 for a non- 
relativistic shock. In Fig. Q we display the 0i dependence of 
equations (6) and (10), in the case of r = 4 for the strong 
shock limit and x = 10 compatible with the assumption of 
weak turbulence (x 3> 1). The upper dotted curve denotes 
the acceleration time without mirror effects normalized by 
that for a parallel shock: 

tA = 4(0i)/4(0i = o°). 

The favoured reduction of the acceleration time for 0i 7^ 0° 
stems from the effects of diffusion anisotropy. It is noted 
that, for 0i 7^ 0°, the reduction of the shock normal compo- 
nent of the particle velocity, v n i, increases the mean res- 
idence time ti, and thus degrades the shock-crossing. As 
seen in Fig.l, this effect dominates the anisotropic diffusion 
effect (coming from the smaller « n i) for 0i > 80°, where 
tA changes to an increase. In the limit of 0i — > 90°, tA 
diverges, because of the related factor (COS01)" 1 in equa- 
tion (10) (see Section 2.3). Hownever, note that the present 
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Figure 1. The acceleration times for an oblique shock normalized 
by that for a parallel shock, as a function of the magnetic field 
inclination 81 (degrees) with respect to the shock normal direc- 
tion. Here, a strong shock (r = 4) and weak turbulence (x = 10) 
have been assumed. The normalized time including the effects of 
both mirror reflection and diffusion anisotropy tA is represented 
by a solid curve and that including solely the effects of diffusion 
anisotropy t\ by a dotted one. The reduction of the time-scale 
indicated by the arrow is ascribed to mirror effects, which are 
pronounced in the large-Si region. 



calculations are physically meaningful only for inclination 
0i ^ cos - (Ui), as mentioned in Section 2.2. For compari- 
son, the normalized acceleration time including the effects 
of both the mirror reflection and diffusion anisotropy, 

t\ = t A (6i)/t A (0i =0°), 

is also plotted (solid curve). That is, the further reduction of 
the acceleration time (from dotted to solid) for 0i 7^ 0° can 
be ascribed to the effects of mirror reflection. It is found that 
the contribution of diffusion anisotropy is larger for smaller 
0i, whereas, for larger 0i, mirror effects play a dominant role 
in reducing the acceleration time. At 0i = 90°, instead of 
the mirroring, the shock drift mechanism actually dominates 
(e.g., Webb, Axford & Terasawa 1983), violating the present 
formalism. 

In order to give a clear presentation of the results, we 
define the contribution rate of the mirror effects as 

Rt = |4 — £a|/4 x 100 (per cent). 

Note that this rate is dependent on 0i, x, and r, and inde- 
pendent of the other parameters. In Fig. [5] for x = 10, we 
plot R t as a function of 0i, given r as a parameter in the 
range 1 < r ^ 4. We mention that, for x 10, the values 
of Rt do not change very much over the whole range of 0i. 
For example, in the case of r = 4, the difference in the R t 
values between x = 10 and 100 is at most 4.3 per cent at 
0i = 79° (not shown in the Figure). In the special case of 
0i = 0° (parallel shock), the effects of both mirror reflection 
and anisotropic diffusion vanish, so that R t — per cent 
irrespective of the compression ratio. As 0i increases, the 
values of R t increase monotonically, and reach nearly 100 
per cent at 0i ~ 90° (quasi-perpendicular shock). As would 
be expected, the contribution of mirror effects is larger in a 
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Figure 2. The contribution rate of mirror effects Rt (per cent) 
against the magnetic field inclination £>i (degrees) for x = 10. The 
cases with the shock compression ratio of r — > 1+ and r = 2, 3 
and 4 are denoted by thin solid, dashed, dotted and thick solid 
curves, respectively. 

stronger shock. For the r = 4 case, the R t value reaches 50 
per cent at 9\ ~ 50° and exceeds 80 per cent at 9\ ~ 74°. 
On the other hand, for r — > 1+ (weak shock limit), R t = 50 
per cent can be achieved at 9\ — 60°. The difference in R t 
for the strong shock case from that for the weak shock case 
is more pronounced in relatively low inclination angles. It 
is also found that, in the range 9i ~ 50° — 70°, the R t val- 
ues merely vary slightly for r 2. Anyway, we can claim 
that the mirror reflection is effective in quasi-perpendicular 
shocks. 



4 CONCLUSIONS 

For a non-relativistic, fast-mode oblique shock accompanied 
by weak MHD turbulence, we have quantitatively revealed 
the contribution of magnetic mirror effects and anisotropic 
diffusion effects to the reduction of the acceleration time 
of cosmic ray particles. We found in particular that, in 
the strong shock limit, for a magnetic field inclination an- 
gle (to the shock normal) of 9i > 50°, mirror effects con- 
tribute dominantly to the reduction of the acceleration time; 
whereas, for 9i < 50°, anisotropic diffusion effects contribute 
dominantly to that time. In the small-#i region, the contri- 
bution rate of mirror effects is found to be small, but sen- 
sitive to the shock strength, such that a larger shock com- 
pression leads to a more enhanced contribution rate. 

While these consequences can be directly referred to 
the study of oblique shock acceleration in space and helio- 
spheric environments, one should be a little more careful for 
an application to other objects, including supernova rem- 
nants. We also remark that the perpendicular diffusion of 
cosmic rays is still not well understood in many astrophysi- 
cal aspects. In a common approach, the diffusion coefficient 
is related to the spectral intensity of magnetic field fluctu- 
ations. For example, spectral analysis of fluctuations in the 
solar corona shows that, in the region of the heliocentric ra- 
dius of 37?0 < R < 6Rq, the power-law indices can be fitted 



by v ~ 1.6, which can be approximated by 5/3 for the Kol- 
mogorov turbulence, and in the 6Rq < R < 12-Rq region, 
v ~ 1.1, whic h can be approximat ed by 1 for the Bohm 
diffusion limit (IChashei et al. 2000). In interplanetary gas 
clouds (extrasolar planetary systems), the power spectrum 
has also been characterized in analogy with the descrip- 
tion of Kolmogorov turbulenc e iWatson. Wiebe fc Crutcherl 
120011 IWiebe k. Watsonll200ll) . Although fluctuations have 
been confirmed in such various objects, this does not always 
mean that the estimated values of x are pertinent to the 
present scheme. 

In young shell-type supernova remnants, especially, 
strong amplification of downstream magnetic field has been 
confirmed by recent observations (e.g. Bamba et al. 2003, 
Vink & Laming 2003). Filamentary structures of hard X- 
rays are interpreted as evidence of shock modification pre- 
dicted by the non-linear theory (e.g. Berezhko, Ksenofontov 
& Volk 2003, Berezhko & Volk 2004, Volk, Berezhko & 
Ksenofontov 2005), and analytic al studies and numerical 
simulations of plasma instability (Bell 2004). The expected 
x values of these sources are of the order of unity, which 
arguably leads to effective acceleration, though the present 
formalism breaks down in the diffusive limit of x = 1 , where 
the magnetic field inclination fluctuates largely, to be am- 
biguously defined. Moreover, it is not appropriate to use the 
unperturbed trajectory of the guiding centre drift motion as 
an approximation while particles are reflected at the shock 
surface by magnetic mirroring. The relevant issues, depart- 
ing from the free crossing approximation, are beyond the 
scope of this paper. 
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